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Abstract 



X 

;h ' We point out that, contrary to some recent claims, there is no intrinsic long-distance un- 

certainty in perturbative calculation of the QCD effects in the tt and tb loops giving the 
electroweak corrections proportional to mf. If these corrections are expressed in terms of 
the "on-shell" mass rrit, the only ambiguity arising is that associated with the definition of 
the "on-shell" mass of a quark. The latter is entirely eliminated if the result is expressed in 
terms of rrit defined at short distances. Applying the Brodsky-Lepage-Mackenzie criterion 
for determining the natural scale for normalization of as, we find that using the "on-shell" 
mass makes this scale numerically small in units of rrit. Specifically, we find that by this 
criterion the first QCD correction to the 0{mf) terms is determined by as^^^{0-l5mt). 
Naturally, a full calculation of three-loop graphs is needed to completely quantify the scale. 



The present remarkable statistics and accuracy of the LEP data at the energy of the 
Z peak calls for a significant theoretical precision of calculation of the electroweak loop 
effects with the goal of sensing the effects of Higgs boson and/or new physics once the top 
quark mass is known (for a recent review see e.g. Ref. [|^). Due to the large mass of the top 
quark, the effects of the tt and ih loops, proportional to m\^ have to be calculated including 
the QCD corrections. Specifically, the leading in the limit mf/m'^^ ^ 1 electroweak 
corrections, related to the W and Z propagators, are universally determined^ by the finite 
difference of the longitudinal parts of the Z boson and the W boson vacuum polarization 
insertions at g^ — > 0. When expressed for definiteness in terms of the correction to the 
electroweak p parameter the effect of the heavy quark loop is given by 
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in the lowest a and the firsts orders in as- The present accuracy of the data already 
makes necessary a quantitative understanding of the magnitude of the as term as well as 
of the higher QCD corrections. 

Recently a doubt was castD B a on the calculability of the higher QCD effects in 
the expression (|l]) in terms of the QCD coupling, normalized at distances of order m^^: 
O-sifnt). The reason for this doubt arises in a calculation of the vacuum polarization 
at q^ ^ 0, i.e. far below the ti and th thresholds, using the dispersion relations, which 
involve integrals over the spectral densities of the physical states, containing ti and th 
near and above the corresponding thresholds. The resonances and the continuum states 
near the threshold are governed by the perturbative and non-perturbative dynamics at 
long distances, which are much larger than m^^, hence the doubt about the calculability 
of the vacuum polarization at g^ -^ in terms of as{fyit). Here we point out that this 
doubt is ungrounded and that the solution of this problem is known long ago: at least 
since the development of the QCD sum rules for charmoniumlii . Moreover, this is exactly 
the central point of the QCD sum rules, that though each individual hadronic state is 
governed by long-distance dynamics, the dispersion integrals over these states, which give 
the vacuum polarization far below the threshold, are determined by the short- distance 
QCD dynamics. This point, which was also emphasized in a recent paper |p, is further 
discussed below in the text. 

The long-distance QCD effects however produce a certain effect on the expression (|1|) 
through the conventions associated with it. Namely, the result in eq.(|ip is written in terms 



of the on-shell mass rrit, which as discussed in this paper effectively lowers the appropriate 
normalization point for as in eq. (|1|) through the contribution of the near-mass-shell region 
to the evolution of the quark mass from the mass shell to distances of order mj~^, which 
are relevant in the loops. In connection with this observation it should be emphasized 
that this effect by no means makes the perturbative calculation uncontrollable: the notion 
of the on-shell quark mass is consistent in any finite order of perturbation theory in QCD 
and the discussed effect manifests itself in a numerical, rather than parametrical reduction 
of the appropriate normalization momentum scale for as in units of rrit. Equivalently, this 
implies that if the QCD corrections to the 0{m^) electroweak corrections are expressed in 
terms of asimt) the coefficients of higher terms should be unnaturally large. Beyond the 
perturbation theory there is the known problem of defining the on-shell quark mass, which 
amounts to an uncertainty of the order of Aqcd, which effectively places the limit on the 
accuracy of definition of the quark mass. This is the well known case for the b quark 
and this particular uncertainty should be the same for the t quark. For the latter the 
situation is somewhat additionally complicated by the large width of the decay t ^ W b. 
Nevertheless, one can, quite probably, get to the accuracy of defining the "on-shell" top 
mass better than ~ 1 GeV. 

One can notice however that the problem of properly defining rrit is somewhat artificial 
for the electroweak corrections at the Z peak. This problem can be completely eliminated 
by expressing those corrections through the top mass, measured at short distances in some 
other measurable quantity, determined by the short distance dynamics. As examples of 
such quantities one can pick the rrit entering the electroweak corrections to the Z —^ bb 
decay rate through the tiW triangle or the total width of the decay t ^ W b, either of 
which is determined by distances of order rraj"^, or any other measurable mt-dependent 
quantity of the same nature. In other words, the relation between the electroweak cor- 
rections of the type in eq.(|l|) and the quantities like the total decay rate of top should 
not contain unnaturally large or small numerical coefficients in units of rrit in the scale of 
normalization of as- 

Naturally, to completely fix the normalization point for as in eq.(|l|) one needs a full 
three- loop calculation in the order a^^, which has not been done yet. However as pointed 
out some time ago by Brodsky, Lepage and Mackenzie^ (BLM) one gets an appropriate 
estimate of the normalization point and thereby makes the coefficient of the higher order 
term reasonable by evaluating the lower order graphs with an explicitly running coupling 
constant. Formally, this corresponds to tagging the dependence of the higher loop term on 



the number of light quark flavors rij and then shifting this dependence in the combination 
60 = 11 — 1^/ into the definition of the normalization point of as in the lower term. Here 
we apply this procedure to the calculation of the electroweak correction in eq.(|I]), and 
find that the BLM criterion gives the normalization point of as in the MS scheme as low 
as 0.15 mt- The appearance of the small coefficient 0.15 is mainly due to the usage of the 
on-shell top quark mass. 



For the practical calculation we use the simple factli-J' that when calculating the 
correction in eq.(0) in the limit m'^/mf <C 1 one can neglect the masses of gauge bosons 
altogether, which is equivalent to setting the electroweak gauge couplings to zero (notice 
that the expression in eq.(||) contains as the overall factor only the top quark Yukawa 
coupling ht). Therefore the quantity of interest can be expressed in terms of the dynamics 
of only the scalar sector coupled to the t and b quarks. In these terms the correction Ap 
is expressed through the vacuum polarization by the quark density operators coupled to 
the Goldstone bosons x° and x^: 



Po(g') = -2z J{0\T{{i{x)^,t{x)){i{x)^,t{x)))\0)e''^^d'x , 

P±iq') = t f{0\T{iiix){l-^,)bix))mx)il + ^,)tix)))\0)e'^^d'x (2) 



(the mass of the b quark is entirely neglected throughout this paper). The electroweak 
correction Ap is found aslii!) 

Ap=^m?(P4(0)-P^(0)) , (3) 

where P'{q'^) = dP{q^)/dq\ 

Consider now a calculation of the difference of the derivatives of the vacuum polariza- 
tion operators in eq.(^ by the conventional Feynman diagram technique. The one-loop 
graphs (no additional gluons) give logarithmic divergence for each of the derivatives, but 
their difference is finite and the result is the leading term in eq.(|l|). Originally the two- 
loop graphs, which lead to the 0{as) correction in eq.(l), were calculated^ by using 
the dispersion relations. Here we concentrate on the direct calculation of these graphs 
by Feynman's technique, which transparently reveals the structure of Euclidean distances 
contributing to the first QCD correction. For this purpose we first consider the inte- 
gration over the quark loop, which allows to represent the result as an integral over the 
gluon momentum. One can notice that starting with the first order QCD correction the 



integration over the quark loop is finite for eacli of the derivative terms in eq. (j^) . Also at 
q"^ ^ there is no obstruction to the Wick rotation. Thus we find in terms of integrals 
over the Euclidean momentum of the gluon the expressions for the first QCD corrections 
to the derivatives of the vacuum polarization operators (0): 



(5«P4(0) = —Jw^l — ] asd{k')k 



rrit J ml 



5«P^(0) = -i^/^o(-)«.rf(fc^)P^ (4) 

47r'^ J \mt mf 



with d{k'^) = l/k"^ being the transverse gluon propagatoiQ and the weight functions w± 
and Wo, determined from the quark loop integration, are given by 

w^ix) = 2 ^yl^^f arctanh [-^y^^j '^-—f^ ln(l + x^)+x^ \nx'-A 

(5) 
and 

Wnix) = 12 — ; ^TTTTTT arctauh , | (6) 

Correspondingly the difference, entering eq.(B|), can be written in the form 



with the resultant weight function 



W{X) = W±[X) — Wo[X 

36 + 70 a;2 + 48 x^ + 12 x^ + x^ 



arctanh 
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x(4 + x2)3/2 \V^+X, 

• ^ + "^'^^ ln(l + x^) + xMn a;2 - 2 ^-^^ (8) 



x^ A + x^ 

At large x the function w{x) has the asymptotic behavior w{x) = 3x~^ + 0{x^^), so 

that the integral in eq.(^ is logarithmically divergent. This divergence is regularized, 

once the gluon propagator is regularized by any standard procedure. Here we do not need 

to specify the regularization procedure, since it can be easily noticed that the expression 

(^) for the measurable quantity Ap contains also the factor mf and the divergence in 

^The expressions (0) are gauge-invariant, hence it is only the transverse gluon propagator, which is 
contributing. 



eq.(l^) is the same as in the renormahzation of m^. More specifically the as correction to 
Ap in terms of the top quark on-shell mass is determined by 

S H (^±(0) - ^o(O))) = m't (- S(m) + 5(P4(0) - P^(0))) , (9) 

where S (7 ■ p) is the one-loop top quark self-energy. The on-shell value of the self-energy 
can be written as an integral over the Euclidean k"^ of the gluon in the loop as 

- S(m) = :^Js(-]as d{k')e^ (10) 



rrit ^ ^ 47r^ J \mt I "^ ^ ' m? 



with the weight function 



,, a;^ + 2x2-8 a;2 
■<'') = ^ 7T=T~^ ■ ^^^) 

Therefore the final expression for the first as correction to Ap (eq.(^) in terms of the 
on-shell mass mt is proportional to the finite integral 

rco 7J-2 

/ {w{x) + s{x)) dx^ = 1 , (12) 

JO 3 

which reproduces the known resultlQ] in eq.(|Tp. 

However, the structure of the integral in eq. (|l2|) deserves a closer look. Indeed, as 
one can see from the plots of the weight functions s{x) and w{x) shown in Fig.l, the 
integral is significantly contributed by the region of small x due to the —2/x behavior of 
the function s{x) in that region, i.e. at k <^ rrit. This behavior is clearly a consequence 
of the fact that the on-shell mass rrit is chosen as the parameter in the electroweak loop. 
As to the weight function of the electroweak loop itself, w{x), as is clearly seen from the 
plot, it is completely dominated by the region x > 1, i.e. k > rrit, and thus it displays 
practically no sensitivity to long distances. 

The significance of the contribution of the region of small Euclidean k can be quantified 
in this calculation by applying the BLM criteriono for the normalization point of as- 
The BLM procedure amounts to replacing the bare gluon propagator asd{k'^) = as/k'^ 
by the one with the running coupling constant: as^{k)/k'^, where as^{k) is the effective 
coupling constant in the potential between two infinitely heavy quarks. The running 
constant as^{k) in the Coulomb gauge corresponds to including the vacuum polarization 
insertions in the propagator of the Coulomb gluon. In particular the BLM procedure 
applied to the calculation in the first order in as correctly reproduces the dependence of 



the coefficient of the next term as^ on the number nj of hght quark flavors, since these 
enter only through the loop insertion in the gluon propagator. However, this procedure 
additionally combines the nj dependence into the factor 60 = 11 — |n/, which is the ffist 
coefficient of the QCD /3-function. In numerous examples'^' LLU applying this criterion to 
the choice of the normalization point for the coupling constant removes large coefficients 
in the subsequent terniQ 

To apply the BLM procedure to the calculation of Ap we use in the integrals in eqs.(|^, 
ToD the one-loop effective coupling: 



as dik') -. 1 «5^(m,) (^1 - ^iMh ln(A:/m,)) . (13) 

Then as times the integral in eq.(|T2p is replaced by 



1 ln(x) {w{x) + s{x)) dx^ 





-as^'im) ( V + 1 1 ( 1 + 1.034 ^^^^gil^) ^ -as^i0.355m,) (^ + l) (14) 

which fixes the normalization point in terms of as^ (the integral with the factor ln(x) was 
calculated numerically.) Furthermore, the normalization of the effective coupling as^ (k) 
is simply relatedi to that of the 05^ as^{k) = as'^ie'^'^k) ^ ^5^^(0.435 A;). 
Therefore from eq.(|T^ we find that within the BLM scheme the effective coupling, en- 
tering the ffist QCD correction in eq.(P is a5*^'^(0.154m(). It is clear, however, that 
to completely quantify the magnitude of the QCD correction to Ap a full three-loop 
calculation of the terms a^^ is needed. 

It is clear that a similar calculation, operating only with Euclidean-space integrals over 
the momenta of the gluons can be in principle performed in higher orders of the QCD 
perturbation theory, thus making it free from the long-distance uncertainties. Exactly 
this point was discussed at length in connection with the QCD sum rules in the papers |^ 



and also in a later review [jT2|. Here we would like to point out a specific loophole in the 

reasoning of the recent papers'^ B El, which state that through the dispersion relations the 

contribution of the near-threshold region to -P'(O) makes the latter quantity sensitive to 

^The only known cases, where large numerical coefficients are not removed by this procedure are those 
associated with the annihilation of heavy quarkonia in QCDo LLiJ, which parallels a similar behavior of 
the three-photon annihilation of ortho-positronium in QED. 



the long-distance dynamics. We also disagree with the argumentation of Ref.|jT3|, where 
it is argued that the threshold effects are numerically small. We insist here that the 
long-distance effects are small parametrically, i.e. suppressed by powers of Kqcd/^^i at 
the non-perturbative level, and are absent altogether in any finite order of perturbation 
theory in QCD. 

The reasoning in those papers is as follows. The difference of the derivatives of the 
vacuum polarization operators in eg.® can be written in the form of the dispersion 
integral 

Pi(0)-Pi(0) = i/^!±W^d. (15) 

where p±{s) and Pq{s) are respectively the spectral densities of the operators (t (1 — 75) h) 
and i a/2 (ij^t), and the integral is running over all values of s, where the spectral densities 
are non-zero. Consider now the region of s near the ti threshold, where the integral is con- 
tributed by the ti resonances and the very beginning of the continuum, strongly distorted 
by long-distance interactions. Within the perturbation theory the exchange of Coulomb 



gluons between the quark and the antiquark with a small velocity v = Jl — Ami/ s makes 
the QCD effects depend on the parameter as/v rather than as- Therefore at v of the 
order of or less than as these effects should be summed up. The summation amounts 
to using the well known solution of the Coulomb problem, and the net effect reduces to 
multiplying the bare spectral density pg by the Coulomb factor 

F. = '""^/'" , , . (16) 

1 — exp(— 47rQ;5'/3v) 

At f ~ 0:5 the spectral density is of order of as and the size of this region of integration 

in eq. ([T5|) is As ~ Ami ^^^ ~ ^(^s^- Therefore the contribution of the 'Coulomb' region 

above the threshold in eq.(^) is of order as^-, which is the same as that of the under- 

the-threshold resonances. The point of the papers 0, H, H is that the as in this effect is 

normalized at long distances: as{m,tv ~ irttas)- Due to favorable numerical factors of tt 

this 0{as^) effect is stated to be a sizeable fraction of the as term in eq.([^). To quantify 

this statement the dispersion integral over the near-threshold region is calculatedo with 

the factor F^. in eq. (p!6D in which the running of the coupling constant is parametrized as 

as{mtv)^. 

■^More precisely, in Ref.|b| a relativistic parametrization is used, which however, does not change the 
main point of the argument. 



The loophole in this argument is that in terms of the running constant as in eq.(|T6D the 
normalization point is not given by nit v- Rather, the proper normalization point is related 
to mtv hy a. function f{as/v): asijntv f{as/v)). This phenomenon is clearly seen in the 
QED calculationli^ of the excitation curve of the r+ r^ at the threshold, including the 
Uehling-Serber running of the QED constant a. Namely, in the region f ~ a the formula, 
obtained by the simple substitution of a by a{'mr v) significantly deviates from the exact 
result. The contribution of the sum over resonances and of the integral over the continuum 
in eq.(|l5|) in higher orders in as contains delicate cancellations!^' li^: starting from order 
as'^ the integral over the continuum partly compensates the contribution of resonance^. 
Moreover, in calculation of the dispersion integral the function Ff. cannot be expanded in 
powers of as-, since this expansion does not converge aX v < as (within such expansion 
the integrals of individual terms would diverge at f = starting from the order (a^/f )^). 
The result of the integration however can be perfectly expanded as a series in as' 111 . 
Thus an approximate parametrization is certainly prone to giving misleading results by 
destroying the correct structure of the spectral density at small v. Therefore we conclude 
that clarifying this point by calculating the spectral density near the threshold is a far 
more complicated problem, than the initial one of calculating the vacuum polarization far 
below the threshold. On the other hand, those detailed calculations of the near-threshold 
region are not needed, since, as discussed before, the 0{m^) electroweak corrections can 
be calculated by the Feynman diagrams entirely in the Euclidean space, in which at no 
step the long-distance uncertainty of the QCD dynamics shows up at the perturbative 
level. In other words, though the integrand in eq.(|l^) is poorly calculable near the ti 
threshold, the entire integral is well calculable within the short-distance QCD. 

As to the non-perturbative QCD effects, these can be understood by adapting the 
results of the discussion of the charmonium sum rules (sections 7.3 - 7.6 of the second 
paper in Ref. |^). The result is that any finite distortion of the quark-antiquark interaction 
at a finite distance tq ^ m~[^ produces only an effect on the ti vacuum polarization at 
q^ ~ 0, which is suppressed by exp(— 2r?7,iro). For instance, one can cut off the Coulomb 
interaction at a radius tq -C {mta)~^ (but still tq ^ ^^)i so that the Coulomb-like 
bound states disappear, and the actual spectral density po would look nothing like that 
determined by the factor F^ in eq. (p!6D . Still, up to the exponentially suppressed terms the 



*The objectionO that these terms in the integral over the continuum can not be negative, since 
the quark and the antiquark are attracting each other, is obviously erroneous: these are small negative 
corrections of order as^ to the positive contribution in the orders as and as'^- 



vacuum polarization at g^ ~ in this situation would be given by the dispersion integral 
in eq.(|T5|) with the perturbative spectral density, i.e. the one containing Coulomb-like 
poles, and the factor Fc above the threshold. Fully appreciating the non-trivial character 
of this phenomenon, we point out that this is a direct consequence of the analyticity of 
quantum amplitudes. The only way, in which the long-distance effects give a contribution 
to the vacuum polarization at g^ ^ is through the 'tail' of the long-distance effects 
at distances of order m^^. In the potential models with a power-like non-perturbative 
potential of the form V{r) = ar"", the effect is proportionallil to arrit (which is the 

action / V{r) dt at distances r ~ mj~^ over the time t ~ m^'^). To evaluate this particular 
effect in QCD there is however no need of invoking model potentials, and the leading effect 
is calculable in terms of the vacuum gluon condensate and its relative magnitude is given 
by {0\7!'asG1^^G'^i^\0)/mf ~ 10"^°. This would be the only non-perturbative contribution 
to the 0{mf) corrections, if these corrections were expressed in terms of the top mass, 
normalized at short distances. However when the corrections are expressed through the 
on-shell mass of top, the relative non-perturbative contribution is that in the on-shell 
mass, i.e. ©(AqcdM*)- 
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Figure Caption 

Figure 1. The weight functions s{x) (left) and w{x) (right) vs. x = k/rrit, in the 
integrals over the Euclidean gluon momentum in the top quark self-energy (eq.(|TUp) and 
in the heavy quark loop (eq. (^ ) . Note the strongly different scale of the vertical axis in 
the plots. 
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